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ABSTRACT 

The recent discovery of two Earth-mass planets in close orbits around an evolved star has raised 
questions as to whether substellar companions can survive encounters with their host stars. We 
consider whether these companions could have been stripped of significant amounts of mass during 
the phase when they orbited through the dense inner envelopes of the giant. We apply the criterion 
derived by Murray et al. for disruption of gravitationally bound objects by ram pressure, to determine 
whether mass loss may have played a role in the histories of these and other recently discovered low- 
mass companions to evolved stars. We find that the brown dwarf and Jovian mass objects circling 
WD 0137-349, SDSS J08205-f0008, and HIP 13044 are most unlikely to have lost significant mass 
during the common envelope phase. However, the Earth-mass planets found around KIC 05807616 
could well be the remnant of one or two Jovian mass planets that lost extensive mass during the 
common envelope phase. 

Subject headings: planetary systems — binaries: close — stars: evolution — stars: low-mass, brown 
dwarfs — stars: individual (KIC 05807616) 



1. INTRODUCTION 
The question of survivability of planets during post- 



some interest over the past years i 


Villaver & Livioll2007l: 


INordhaus et al.ll2010l:ISDieeel2012 


) . Recently, many sub- 



stellar companions have been discovered in close orbits 
around evolved stars (Table [T]). Maxted et al. (2006) de- 
tected a 0.053 Mq brown dwarf orbiting the 0.39 Mq 
white dwarf WD 0137-349 with a 1.9-hour period. A 
likely brown dwarf was discovered in a 2.3-hour or- 
bit around the s u bdwarf B star SPSS J08205-t-0008 
(|Geier et all 120111) . iSetiawan et~all (|2010[ ) discovered a 
Jupiter-mass object orbiting the red horizontal branch 
star HIP 1304 4 with a 16.2-day period. Finally, 
iCharpinet et al.l (|2011l ) reported the detection of two 
nearly Earth-sized planets orbiting the subdwarf B star 
KIC 05807616 with orbital periods of 5.8 and 8.2 hours. 
The existence of these systems for which the radius of the 
precursor of the primary was larger than today's orbital 
separation, suggests that they have gone through a com- 
mon envelope (CE) interaction. During this phase, the 
companion is engulfed by i ts host star's en velope when 
the latter was a giant star (lPaczvns"killl976[ ). 

Hydrodynamics simulations of this evolutionary phase 
have been p erformed with d ifferent numerical techniques 
(ISandauist et al . 1998: P assv et al.|[2012,: Ricker fc TaamI 
I2012f) but questions persist. It is still unclear how the 
envelope of the donor star gets ejected, and why the final 
separations obtained by simulat ions are larger tha n those 
observed in post-CE systems (jPassv et al.l 120121 ). The 



resolution reached in these numerical models currently 
does not suffice to study the complete evolution for a 
substellar companion that may spiral down to very close 
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to the giant's core. 

An important question is therefore whether substellar 
companions can survive engulfment by the giant's enve- 
lope without being totally disrupted. A system with a 
very low mass companion has a very small orbital energy 
budget, so we may also plausibly conjecture that a CE 
in such a system would more likely fail and result in a 
merger, thus destroying the companion. Is it therefore 
possible that substellar companions started off as more 
massive objects, which were then partially disrupted dur- 
ing their engulfment to become the obje cts we see today? 
This s cenario was actually suggested by ICharpinet et al.l 
(|2011f ) to explain the existence of planets around KIC 
05 807616. 

. Murrav et al.l (|199 3') studied the disruption of gravita- 
tionally bound objects by direct action of ram pressure 
in the context of star formation and galaxy evolution. 
However the physics of disruption is scale free, so their 
results also apply to other astrophysical questions such as 
the gas stripping occurring in clustered galaxies (i Nulsenl 
[1982; Mori & Burkcrt 2000) or here, the disruption of 
hydrostatic companions captured by the envelope of a 
giant star. In this contribution we use their theory to 
determine whether, in the observed systems, the substel- 
lar companions could have lost a significant amount of 
mass during the in-spiral phase. We describe the formal- 
ism used in Section[2]and present the results in Section|3l 
We provide a summary in Section 3) 

2. ANALYSIS 

iMurrav et al.l (|1993l ) calculated the conditions for a 
self-gravitating, hydrostatic, isothermal sphere moving 
in an ambient medium to be disrupted by ram pressure. 
The disruption of an object by ram pressure can most 
simply be understood by the pancake model (Zahnl^ 
1992; iKlein et al.lll994| ). in which transverse motions are 
driven by the increased pressure on the front surface o f 
the object compared to the sides. iField fc Ferraral (|19950 
showed that the actual disruption by Kelvin-Helmholtz 
and Rayleigh- Taylor instabilities with wavelengths com- 



2 



parable to the size of the object reproduces the scal - 
ings yielded by this simple model. iMurrav et al.l (|1993[ ) 
confirmed numerically that when self-gravity is strong 
enough to stabilize the Kelvin-Helmholtz instability, the 
sphere is stable against catastrophic disruption. This oc- 
curs when the gravitational acceleration 
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(1) 



where pi, R2 and U are the background density, and 
the density, radius and relative velocity of the compan- 
ion, respectively. We assume the tidal forces exerted by 
the giant star on the companion to be negligible com- 
pared with the self-gravity of the companion - we discuss 
this point in Section 13.4.11 Assuming D = p2 /pi 3> 1 
and taking the gravitational acceleration g w GM2/R2, 
one obtains a critical mass under which disruption occurs 
(their equation 2.8): 
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If the configuration is unstable (M2 < Merit), the com- 
panion is disrupted in a crossing time as the characteris- 
tic timescale for destruction is equivalent to its internal 
dynamical timescale: 



1.4 



R2 



0.1 R, 



/lOO kms"^ 

V u 



1/2 

50 



We approximate the relative velocity of the companion 
by the maximum orbital velocity reached during the in- 
spiral phase, i.e. the final orbital velocity: 
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where Mc ^ M2 is the mass of the primary's core, and a 
is the orbital separation at the end of the CE phase. We 
substitute this relation in Equation ([5]), and obtain the 
critical background density above which the companion 
is unstable: 
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One should emphasize that if the companion is unstable, 
the dominant wavelength acting to destroy it is compa- 
rable to the sphere radius. The destruction is therefore 
global and not quasi-static. In other words, the com- 
panion is destroyed catastrophically, not ablated by in- 
cremental mass-loss due to short wavelength surface in- 
stabilities. Thus stability cannot be reached again once 
disruption starts. 

If the companion is stable (M2 > Merit), such shorter 
wavelength instabilities can still occur and strip mass in 
a more ablative manner, giving a quasi-static mass loss 



rate. We compute the largest unstable wavelength from 
Equations ([T]) and ([4]), again assuming p2 ^ pi, and 
substituting Amax for R2. We find 
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IMurrav et al.l (|1993t ) suggests that the corresponding 
mass loss rate can then be approximated to that oc- 
curring in laminar viscous fiows as described by their 
Equation 3.5: 
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where we have assumed the relative velocity to be of 
the order of the sound speed of the background. Equa- 
tion (O should be used with caution as it is valid only for 
Reynolds numbers Re < 30 ()Nulsen|[l982[ ). In our case, 
the Reynolds number is Re = 2R2U/1', where ly is the 
molecular viscosity. Viscosity in stellar interiors is far 
too small to allow such low Reynolds numbers on these 
length scales. 

However, such low effective Reynolds numbers might 
be reached through the turbulent viscosity. Since the 
largest possible size of the eddies is Amax and their largest 
velocity is the sound speed of the ambient mediu m c.q, 
one ca n use a formalism like that of iShakura fc SunvaevI 
(jl973| ) and parametrize the turbulent viscosity as: 



(8) 



where a < 1 is a free parameter. The effective Reynolds 
number is thus: 
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where we have still assumed that [/ ~ eg. Unfortunately, 
no accepted value of a has been derived for our spe- 
cific case. We can only calculate the minimum turbulent 
Reynolds number -Ret.min = Ret{a = 1) and see how this 
number compares to the range required to apply Equa- 
tion dZl). 

Another approach would be to calculate how much 
work from the drag acts on the companion during its in- 
spiral, and compare it to the companion's gravitational 
binding energy E'bin ^ GM2/R2- Assuming for the cross 
section S2 = 7ri?2 and a drag coefficient Cd, the drag 
force is Fd — \CdPiU^ 82- The smaller the compan- 
ion, the longer the in-spiral phase lasts: for example, the 
0.9 Mq and 0.1 companions comple te about 20 and 
40 orbits during this pha se, respectively (|De Marco et al.l 
[Ml IPassv etaII[20Tl) . We therefore assume that our 
companions complete at least 50 orbits, and slice the gi- 
ant's envelope down to the currently observed separation 
of the system, in iVorbits = 50 equally-spaced radii Ri cor- 
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responding to an average density pi. Assuming C^i = 1, 
the work provided by drag forces is thus: 



(10) 



where Ui = \J GMi/Ri and Mi is the mass enclosed 
within Ri. We will see in Section 13.41 that the calcu- 
lation of Wdrag does not depend strongly on the number 
of assumed orbits. If W^drag ^ Ehm we conclude that 
the companion could have been destroyed if the drag en- 
ergy can be coupled to the interior of the companion 
efficiently, even if it is stable according to Equation ([T]). 

3. RESULTS 

In this section we study each system reported in Ta- 
ble [T] separately and investigate whether the substellar 
components are less massive today than they were prior 
to the CE evolution. 

3.1. WD 0137-349 
Using the MESA stellar evolution code (jPaxton et al.l 



|2011|) . we compute a model for the 0.053 Mq brown 
dwarf companion to this white dwarf and obtain a ra- 
dius i?2 = 0.079 R© (Figure [1]). The mean density of the 
brown dwarf is p2 = 150 g cm~'^ which, substituted in 
Equation ([5]), yields a critical density pi^crit ~ 27 g cm~^. 




Fig. 1 
(dashed 
models. 



0.12 



— Density profiles for the 0.053 M0 (solid black), 10 Mj 
red), 1 Mj (dash-dotted blue) and 0.1 Mj (dotted green) 



We compare this critical density to the density profile 
of the progenitor and see whether the companion will 
plunge deep enough into the primary's envelope to en- 
counter a density high enough to destroy it. The most 
likely progenitor of the white dwarf is a nearly 1 Mq 
main sequence star. We evolve such a model with MESA 
up to the red giant branch until the core mass reaches 
the observed white dwarf mass. At that time, the radius 
of the star was 100 Kq and its total mass was 0.89 Mq 
due to mass loss. The density profile of this red giant 



star is displayed in Figured The current orbital separa- 
tion of the system is a = 0.65 Rq which means that the 
companion has reached a layer of the primary where the 
density pi,max 2.1 x 10"*^ g cm~^. 

At this point, we should emphasize that we have cal- 
culated the maximum density encountered by the com- 
panion using the initial profile of the giant primary. Hy- 
drodynamical simulations suggest that in the deep inte- 
rior, density at a given coordinate decreases with time 
(jPassv et al.ll2012L Figure 12) although the current reso- 
lution does not allow them to give a definitive answer at 
such small radii. Therefore, pi,max should be considered 
an upper limit. 

Pi, max is four ordcrs of magnitude smaller than the 
critical density required for instability so we conclude 
that the companion is stable. Moreover, the mass loss 
rate due to viscous stripping ~ 3.8 x 10~^ Mj yr^^ only. 
Since the dynamical phase of a CE lasts only about one 
year (|Sandauist et al.lll998l: iPassv et al.ll2012[ ). the com- 
panion has not lost any significant amount of mass during 
its in-spiral phase. 




Fig. 2. — Density profile for the 0.89 Mq red giant star progenitor 
of WD 0137-349 (soUd black) and a 0.76 Mq star at the tip of the 
red giant branch (dotted red). 



3.2. SDSS J08205+0008 

SDSS J08205-f0008 is very similar to WD 0137-349 
so we also conclude that the brown dwarf - whose es- 
timated mass ranges from 0.045 to 0.068 Mq - has not 
been affected during the plunge-in phase. 

3.3. HIP 130U 
The minimum mass of th e planet orbiting HIP 13044 
is estimated to be 1.25 Mj (jSetiawan et al.ll2010D . This 
system is intriguing, as a plausible formation channel 
for this system has yet to be found. Indeed, the or- 
bital separation (~ 25 Rq) is much larger than the 
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separations obtained from numer ical simulations or ob- 
servations of pos t-CE systems ()De Marco et al.l 120111 : 
iPassv et"al1l2012D . Also, successful CE interactions usu- 
ally lead to the formation of blue ho rizontal branch s tars 
with very thi n hydrogen-envelopes (iHan et al1l2002| ). 

iBear et al.l (2011) therefore suggested that the he- 
lium flash occurring in a 0.9 M0 metal-poor red giant 
star could ignite hydrogen burning in a thick region at 
the bottom of the envelope. This deposition of energy 
would eventually lead to a dramatic inflation lasting 
about 100 yr, during which the stellar radius could reach 
700 R© . The planet would then enter a "half- failed, half- 
successful" CE during which the orbital separation would 
shrink down to approximately 25 R©, only 0.1 M© would 
be lost by the system and the primary would evolve to 
the red part of the horizontal branch. 

This theory has several problems. First, there is no 
compelling evidence suggesting that red giant branch 
stars might suffer such a large and long-lasting expan- 
sion after the helium core flash. On the contrary simu- 
lations performed with stellar evolution codes such as 
MESA, show that red giant stars shrink after helium 
has been ignited explosively. Secondly, the duration of 
this suggested post-helium-flash expansion is larger than 
the duration of a CE phase. Therefore, the in-spiral of 
the companion could not have been stopped because of 
the primary contracting below its Roche lobe on a short 
timescale. Thus, the entire envelope mus t have been 
lifted up, mostly in the e quatorial plane (jPassy et al.l 
[20TI iRicker fc TaamllMl . which raises thus the ques- 
tion of how a red horizontal branch star could eventually 
form. Finally, the planetary companion would strongly 
interact with falling back material through gravitational 
drag. In that regard, the large orbital separation ob- 
served is still perplexing. For all t hese reasons we b e- 
lieve that the scenario suggested bv IBear et al.l (|2011l ) is 
unlikely. 

A careful investigation of the progenitor of HIP 13044 
seems quite complex and beyond the purposes of this pa- 
per. Consequently we simply assume that the companion 
went through a CE interaction with an unknown progen- 
itor. For the latter, we will consider the same profile as 
in Section 13.1! in order to give an indicative answer on 
the potential survival of the companion. 

For the companion, we consider a 1 Mj model (Fig- 
ure [1]), with a radius R2 = 0.101 R© and a mean den- 
sity p2 = 1.31 g cm^"^. We find the critical density 
for stability pi^crit ^ 6.2 x 10~^ g cm~^. The den- 
sity in the progenitor's profile corresponding to the ob- 
served separations between the subdwarf and the planet 
Pi, max ~ 10~^ g cm~^, almost 4 orders of magnitude 
smaller than pi^crit- We thus conclude that the compan- 
ion is stable. The viscous mass loss rate is also negligible 
(Mv ^ 1.4 X 10~* Mj yr~^) so no mass will be ablated 
from the companion during the plunge-in phase. 

3.4. KIC 05807616 

The most intriguing system is certainly KIC 05807616, 
which has two nearly Earth-mass planets orbiting a sub- 
dwarf B star with an orbital separation of about 1.5 R©. 
In order to constrain the mass of these planets prior to 
their engulfment, we use the same method as previously, 
with four different models ranging from 10 Mj down to 



1 M© (Figured]). For the Earth-sized mo del, we use the 
accepted values for the Earth (see, e.g., iCoxl |2000|) as 
listed in Table [2j We again assume the primary used in 
Section lOI for which we obtain pi,max ~ 10~^ g cm~'^. 

The results are summarized in Table [2] One can see 
that the critical density required to catastrophically de- 
stroy an Earth-sized companion is about as large as the 
maximum density pi,max- In this regard, the 1 M© com- 
panion could still be stable. However Wdrag ~ lO'*^ ergs, 
which is more than an order of magnitude larger than 
the gravitational binding energy of the companion (Ta- 
ble [2]). Wdrag does not strongly depend on the number of 
assumed orbits since changing this number to 20 or 100 
changes Wdrag by a factor 2-3 only. Since the work pro- 
vided by drag forces exceeds the binding energy of the 
companion, we conclude that the Earth-sized compan- 
ion may not survive the engulfment, if the energy can be 
effectively coupled to the planetary material. 

From the same considerations, a 10 Mj companion is 
too massive to be affected in any way during its in-spiral 
phase (Table [2]). For a 1 Mj companion, pi^max exceeds 
Pi,crit by almost an order of magnitude, so the secondary 
could be stable. However, E'bin exceeds Wdrag by only a 
factor of 3, which suggests that a significant amount of 
mass could be removed from the surface of the compan- 
ion. A 0.1 Mj companion has a critical density an order 
of magnitude smaller than pi^rnax- Moreover, Wdrag is 
about twenty times larger than the binding energy of the 
companion. Together, these suggest that an 0.1 Mj ob- 
ject probably would be destroyed. Using a more evolved 
model with a higher core mass Mc = 0.47 M© and 
A'^orbits = 100 (the stellar radius ^ 175 R©) yields sim- 
ilar results (Figure m Table [2]) We thus conclude that 
the progenitor of the planets likely had a mass of a few 
M.T, s upporting the suggestion made b\ lCharpinet et all 

mm . 

3.4.1. Discussion 

A legitimate question is whether mass accretion by the 
substellar compa nion is relev ant. From a Bondi-Hoyle 
analysis (see, e.g.. lEdgail2004 equation 31) the accretion 
rate onto the 1 Mj object is only 

Mace - ""^'^^r^^ = 3.7 X 10-2 Mj/yr. (11) 

We therefore conclude that mass accret ion is negligible 
on the timescale of the CE interaction. IRicker fc TaamI 
(|2008f) also reached this conclusion in the case of more 
massive secondaries. They further argued that the mass 
accretion rate calculated from their simulations is even 
smaller than the one given by Equation ()11|) . mainly be- 
c ause tidal effec t s dora inate the structure of the flow. 

IBear fc Soke^ (I201I) propose a slightly different sce- 
nario for the formation of the two planets around KIC 
05807616. Instead of considering the dynamical effects 
occurring during the CE phase, they suggest that a mas- 
sive giant planet was tidally destroyed after the end of the 
in-spiral phase and that the two planets are remnants of 
the disrupted metallic core. Usi ng their equation 1 wit h 
C'tido = 1, the values given bv iCharpinet et al.l (|201lD . 
and our planetary models (Table [2]), the tidal radius at 
which destruction occurs is Rt = 0.39 and 0.80 R© for 
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TABLE 1 

Orbital parameters. 



Name 



Mc (Mq) 



M2 



a (Rq) 



References 



WD 0137-349 
SDSS J08205+0008 
HIP 13044 
KIC 05807616 
KIC 05807616 



0.39 
0.25-0.47 
0.8 
0.496 
0.496 



55.6 Mj 
47.1 - 71.2 Mj 
>1.25 Mj 
0.440 M0 
0.655 M® 



1.93 hours 0.65 Maxted et al. (2006') 

2.30 hours 0.60 - 0.72 Geier et al. (20TT) 

16.2 days 24.95 Setiawan et al. (2010) 

5.7625 hours 1.290 Van Grootel et al. (2010); Charpinet et al. (20n) 

8.2293 hours 1.636 Van Grootel et al. (2010); Charpinet et al. (20^ 



Note. — Reported are the stellar mass (-A/c), the companion mass {M2), the orbital period (P) and the orbital separation (a) 



the 10 and 1 Mj objects, respectively. These quahtative 
values are still sraaller than the observed separations of 
the planets (see Table [T]) so it is difBcult to conclude 
whether these objects would be tidally destroyed. They 
also suggest that evaporation of the companion due to 
heating from the envelope might be important. They es- 
timate the mass of the planet must be > 5 Mj in order to 
survive evaporation. The evaporation timescale depends 
on the efficiency of heat transfer, which is sensitive to 
mixing between the planetary atmosphere and the stel- 
lar convective envelope. We cannot at present calculate 
this timescale with sufficient precision. We therefore do 
not know whether evaporation could dominate the dy- 
namical effects we discuss in this study. Nevertheless, 
the Eu ass of the progenitor suggested by iBear fc Sokeii 
()2012t ) is consistent with ours. 

4. SUMMARY 

We have considered evolved stars found to have brown 
dwarfs or planets around them in close orbits, suggesting 
that these systems have gone thro ugh a CE inte r action . 
Using the criterion developed by iMurrav et al.l ()1993[ ). 
we have investigated the stability of these substellar com- 



panions against mass loss induced by Kelvin-Helmholtz 
and Rayleigh- Taylor instabilities. 

We have found that the substellar companions ob- 
served in WD 0137-349, SDSS J08205+0008 and HIP 
13044 are very unlikely to have been affected during the 
in-spiral inside their giant progenitors. Therefore, the 
masses of these objects currently observed are likely the 
same as prior to the CE phase. The two planets detected 
orbiting KIC 05807616, however, are not massive enough 
to have survived the engulfment unscathed. We have es- 
timated the mass of the progenitor of these planets to be 
of order 1 Mj. The question that still has to be answered 
is how such low mass companions could have successfully 
unbound the entire envelope of the primary. 
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TABLE 2 

Parameters of the different companion models investigated for KIC 05807616. 



Ma 


R2 (Rq) 


p2 (g.cm 3) 


Pl.crit (g-cm ^) 


-^max / R2 


^^t,min 


Mv (Mjyr-l) 


M^drag (ergs) 


Ebin (ergs) 


10 Mj 


0.106 


11.4 


0.5 


2.0(-3) 


9.9(2) 


0.6 


1.2(43) 


3.3(45) 


1 Mj 


0.101 


1.31 


6.0(-3) 


1.7(-1) 


12 


45 


1.1(4,3) 


3.4(43) 


0.1 Mj 


0.080 


0.27 


1.5(-4) 








7.0(42) 


4.3(41) 


1 Me 


9.16(-3) 


5.52 


8.8(-4) 








1.1(41) 


3.7(39) 


10 Mj 


0.106 


11.4 


0.5 


1.0(-3) 


2.0(3) 


0.15 


5.3(42) 


3.3(45) 


1 Mj 


0.101 


1.31 


6.0(-3) 


8.4(-2) 


24 


12 


4.8(42) 


3.4(43) 


0.1 Mj 


0.080 


0.27 


1.5(-4) 








3.0(42) 


4.3(41) 


1 M® 


9.16(-3) 


5.52 


8.8(-4) 








4.9(40) 


3.7(39) 



Note. — Values obtained with the 0.89 Mq (top) and the 0.76 Mq (bottom) models. Quantities relative to viscous stripping 
are not calculated when the object is found unstable. 



